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ABSTRACT 

Context. Low mass X-ray binaries were found to have very low frequency breaks in their power density spectra below 
which the power density spectra are nearly in white noise structure and at higher frequencies they approximately follow 



the Py (XV 



law. 



Aims. To determine whether high mass X-ray binary power density spectra have similar properties and the findings for 
low mass X-ray binaries are also valid for high mass binaries, we analyzed the time series of high mass X-ray binary 
sources produced by All Sky Monitor of Rossi X-ray Timing Explorer. 

Method. We obtained the power density spectra of the high mass X-ray binaries using the cosine transform of autocor- 
relation function. 

Results. We identified break frequencies for seven sources, namely OAO 1657-415, SS 433, Vela X-1, SMC X-1, 4U 1700- 
377, GX 301-2, and LMC X-1. The normalized break frequencies with respect to the orbital frequency {f break / fortu) 
for sources OAO 1657-415, SS 433, SMC X-1 and LMC X-1 are consistent with those of Roche lobe overflow systems. 
The other high mass X-ray binary systems. Vela X-1, GX 301-2, and 4U 1700-377, however, have larger break frequency 
ratios, fbreak/ forb, which are indicative of short viscous times. These are all wind-accreting sources and the stellar winds 
in the systems allow the formation of only short radius discs. Consequently, we qualitatively distinguished the Roche 
lobe overflow binaries from the wind accreting system by comparing their normalized break frequencies. 
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1. INTRODUCTION 



X-ray binaries have X-ray flux variations on a broad range 
of timescales. In a binary system tlie transfer of matter 
through an accretion disc is considered to be dominated by 
viscous forces. These modify the disc structure on a viscous 
timescale. Hence, the response of the disc to any external 
change occurs on this longest timescale of the disc. If the 
frequency of a variation is lower than tmsc{Rd)~^ ~ fbreak, 
the corresponding X-ray flux variations are considered to 
be uncorrelated; therefore, the corresponding power density 
spectrum is expec ted to exhibit a white no ise behavior at 
these frequencies (jGilfanov fc Arefievll2005[ ). 



iGilfanov fc ArefievI (|2005D found that the power density 
spectra of some persistent low mass X-ray binary systems 
have characteristic shape: white noise up to a certain fre- 
quency value after which a power law of P (x appears 
where a ^ 1.3. 



In this study, we find significant break frequencies in 
OAO 1657-415, SS 433, Vela X-1, SMC X-1, 4U 1700-377, 
GX 301-2, and LMC X-1. In the following section, we ex- 
plain the selection criteria and our sample of X-ray binaries. 
We also the describe the techniques we used to obtain the 
power density spectra, and then determine the break fre- 
quencies together with the power law indices. Discussion of 
our results is given in section 3. 



2. DATA 8i ANALYSIS 

2.1. Data 

We used data gathered by the All Sky Monitor (ASM) in- 
strument of Rossi X-ray Timing Explorer (RXTE) obser- 
vatory, covering the period from 1996 to 2010. The ASM 
instrument operates in the 2-12 keV energy range, per- 
forms flux measurements for over ~ 500 X-ray sources, and 
scans 80% of the sky every 90 minutes. Each flux measure- 
ment has a duration of ~ 90 seconds. Some gaps some- 
times appear in time as long as a few months in the light 
curves of ASM sources because of navigational constrains 
and the appearance of very bright transient sources. The 
dwell-by-dwell light curves were retrieved from the public 
RXTE/ ASM archive at HE AS ARC. 

We derived the power density spectra (PDS) of the 
sources from the cosine transform of the discrete autocorre- 
lation function. This technique has advantages if the tirn e 
series are unevenly sampled ^Gilfanov & ArefievI |20Q5|). 
We also estimated the PDS using HEASARC software 
" pow spec" with Miyamoto normalization ([Miyamoto et all 
Il992[ ) in our data analysis. Both estimations of PDS give 
consistent results. 



2.2. Sources 

We constructed our sample after evaluating the PDS of light 
curves for 45 high mass X-ray binaries with known orbital 
periods. 
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Fig. 1: Long-term light curves of the high-mass X-ray binaries 
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Table 1: The binary parameters for LMXBs from the ASM sample 
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Source 


Mi(Mq) 


M2{Mq) 


Porb(days) 


q 


D(kpc) 


Ref. 


UAUiDO (-4io 


1 A 


1 A 1 Q 

14-18 


1 n /I /I 'id 
iU.44oD 


iU-iz.y 


D.4 zt i.Z 


i 


SS433 


~ 9 


~ 30 


13.082 


- 3.33 


5 


2 


Vela X-1 


2.27 ±0.17 
1.88 ±0.13 


27.9 ± 1.3 
23.1 ±0.2 


8.964368 


12.29 ±0.35 
12.29 ±0.74 


1.9 ±0.2 3 


4 


GX 301-2 


1.85 ±0.6 


39-53 


41.5 


21.74 ± 6.62 


3-4 


5 


SMC X-1 


1.06 ±0.1 


17.2 ±0.6 


3.8923 


16.23 ± 2.10 


61 


6 


LMC X-1 


10.91 ± 1.41 


31.79 ±3.48 


3.90917 


2.91 ±0.70 


~ 50 


7 


4U 1700-377 


2.44 ± 0.27 


58± 11 


3.412 


23.77 ±7.14 


1.8 


8 



Notes. Ml and M2 are masses of the primary and secondary companions, respectively. Port is the orbital period, a is the binary 
separation, q is the mass ratio M1/M2, and D is the distance of the system. The mass of the pulsar of OAO 1657-415 has 
not been determined to be a certain value and we assumed it to be 1.4M0 in our calculations. The two sets of parameters for 
HMXRB corresponds to incli nation angles i = 70.1 deg a nd 9 0deg, respective l y. Re f. re fers to References for th e p arameters: (1) 
iChakrabartv et all (Il993l): (2) iCherepashchuk et all (120051): (3 ) lSadakane eT^ (jigsg): (4) iQuaintrell et al.l 12003 ): (5) lKaper et all 
l|2006t ): f6) lvan der Meer et all l|2007D : f7) lOrosz et al.l (|2009l ): ('8) IClark etalf |2002 ) 



(i) We selected only persistent sources that are in the 
same state during the observation period. Twenty- 
five HMXRBs in our sample are Be transients that 
are not persistent X-ray sources and display outbursts 
because of the extremely wide eccentric orbit of neu- 
tron stars around their companions. Owing to the 
short outburst time span relative to the orbital pe- 
riod, we do not see breaks in their PDS. Outburst 
phases are the cause of substantial red noise in their 
PDS. Two sources Cen X-3 and Cyg X-1 have low fre- 
quency breaks in their PDS. However, these sources 
have significant signatures of high-low states in their 
light curves. These state transitions contaminate the 
long-term power density spectra, so that fake breaks 
may appear. To see these artificial effects, we carried 
out simulations as described below in (iii). In addition, 
we refer to Fig.[2]for the comparison between the PDS 
obtained from the real light curve and the simulated 
light curve. While creating the fake light curve, we 
took into account the state transitions of the source. 
Although we produced this light curve from a PDS 
obeying a power law trend (i.e. no breaks), it is clear 
from the figure that the simulated PDS resembles the 
real PDS. Hence, we estimate that the breaks in these 
sources are associated with the occurrence of high- 
low intensity states rather than the viscous time of 
the disc. 

(ii) The sources with very low count rates, < 1 count/sec, 
such as LMC X-4 or 4U1538-52, have white-noise- 
dominated power density spectra at low frequencies. 

(iii) In our analysis, we found that the remaining seven 
sources had low frequency breaks in their PDS to- 
gether with powe r-law indices ~ 1.3 c onsist ent with 
the LMXRBs of iGilfanov fc ArefievI (pOOSl ). These 
are OAO 1657-415, SS 433, Vela X-1, SMC X-1, 4U 
1700-377, GX 301-2, and LMC X-1. These values 
for a in dices imply the existence of accretion discs 
(|Lvubars kii 1997). In Figs. [Hand Hand Table [H light 
curves, PDS and orbital parameters of the sources are 
presented, respectively. We marked the orbital and 
super-orbital periods of the sources on the power den- 
sity spectra. Although these modulations are quite 
clear as peaks at the corresponding frequencies, it 
is not possible to determine their values from the 



PDS accurately since we use geometrical binning and 
the resolutions of the PDS are low. Table 2 contains 
the break frequencies and power-law indices of the 
sources. 

Power spectra of light curves of the sources show that 
there are periodic modulations in the luminosity of the 
sources that are caused by both the orbital motion 
of the system and the quasi-periodic/periodic mod- 
ulation due to the precession of accretion disc. To 
see whether the low frequency breaks, which we saw 
in the PDS of the HMXBs, are the results of such 
periodic/quasi-periodic variability in the X-ray flux, 
we carried out simulations using the orbital period 
and precession periods of each source. In simulating 
the light curves of each source, we took Fourier ampli- 
tudes from the real power spectra, and for the phases 
of Fourier amplitudes, we used Gaussian white noise 
with unit variance. We sampled the light curves of 
each source as if observed by ASM such that their 
power densities are consistent with the power law 
P ^ f~^'^ ■ We then convolved these time series by in- 
troducing periodic signals for both the orbital motion 
and precession period if the source had a periodic- 
ity. Consequently, using the simulated light curves we 
constructed the power density spectrum. In our simu- 
lations, we did not see any breaks due to the orbital or 
precession modulation of light curves. The simulation 
result for Vela X-1 is presented in Fig (21 
Other short-term variations, such as flickering and 
magnetic field effects on the inner disc, appear at high 
frequencies, i.e. higher than our Nyquist frequency 
> 1.0 X 10~* Hz, so they do not 'pollute' the low 
frequency power density spectra we analyze. These 
short-term variations are the subje ct of independent 
resear ch. (For magnetic effects, see iRevnivtsev et"aLl 
(120091) .) Therefore, aU of the seven sources were cho- 
sen for further analysis in this work. 



The primary goals of this work are to study the low- 
resolution PDS of HMXBs and both resolve low frequency 
breaks and distinguish them on the basis of their accretion 
properties. According to the standard theory of the alpha 



3. DISCUSSION 
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Fig. 2: On the left side, the simulation results of Cen X-3 is presented. The top figure is the power density spectrum in power law 
with orbital modulation, the middle one is the same as the former but the state transitions are taken into account, and the bottom 
figure is the real PDS of Cen X-3. On the right side, the simulation results of Vela X-1 is presented. The top figure is the power 
spectrum only in power law, the middle one is the PDS in power law with orbital modulation, and the bottom figure is the real 
PDS of Vela X-1. It is clear from the figures that sinusoidal orbital modulation does not affect the trend of the power law, but the 
state transitions cause an artificial break in the PDS. 



disc model (Gilfanov & Arefiev 2005, Shakura & 
Sunayev 1973), the viscous timescale of accretion disc (in- 
verse break frequency) Uisc = '^1 fvisc ^ hreak is pro- 
portional to the orbital period and this proportionality is 
given by 



27ra 



Rd 



Rd 

a 



-3/2 



(1) 



where a is the viscosity parameter and H/ R is the ratio 
of disc thickness to outer disc radius, q is the mass ratio 
M1/M2, and a is the binary separation. 

In Fig. m we presented the correlation between the 
break frequency and orb ital frequency of HMXB s as well 
as the LMXBs observed (|Gilfanov fc Arefiev|[2005l) for var- 
ious H/R and q values, respectively. 

In Fig.[5]-a, we plotted the normalized break frequencies, 
fbreak/ forb, against the orbital periods. In this figure, it can 
clearly be seen that, for Vela X-1, 4U 1700-377, and GX 301- 
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OAO 1657-415 









Fig. 3: Power density spectra of high mass X-ray binaries with the best power-law models (The PDS of Vela X-1 is presented in 
Fig[2]at the bottom right corner) 



2, the normalized break frequency is greater than for other 
sources (see Table 2). For this group of sources, normalized 
break frequencies are in the range 0.76-2.56, whereas for SS 
433, SMC X-1, GX 301-2, and OAO 1657-415 these values 
are in the range 0.076-0.16. The ratios of the normalized 
frequencies of the HMXBs in the former group to those in 
the latter group range between 5 and 33. If the normalized 
break frequencies are associated with viscosity times, then 
Vela X-1, 4U 1700-377, and GX 301-2 should have tiny discs 
with respect to the other sources. 

There should be some reasons for this distinction re- 
garding the properties of accretion processes in these bi- 
naries. Since we are investigate high mass X-ray binaries, 
there are two main accretion mechanisms: Roche lobe over- 



flow and stellar wind accretion. Roche lobe overflow leads 
to persistent disc formation around the compact object. 
On the other hand, thin accretion discs may appear in the 
course of wind accretion. Among our sources, SS 433 and 
SMC X-1 have been determined to be Roche- lobe overflow 
binaries. Vela X-1, GX 301-2 and 4U 1700-377 involve wind- 
driven accretion mechanism, and for the rest, namely OAO 
1657-415 and LMC X-1, there are remaining uncertainties 
about the nature of accretion mechanisms since they have 
been observed to show signs of both Roche-lobe and wind 
accretion. To understand these categories more clearly, we 
investigate in the following section the accretion properties 
of each source. 
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Table 2: The parameters of power spectrum approximation 



Source 


(count / sec) 


Lx 

(lO^^ergs/s) 


/i)reafe(Hz) 


Q 


f range (Hz) 


f break/ forb 



OA01657-415 


2.25 


0.3 


-| 74-1-0.30 

^■'^-0.30 


X 10" 


-7 


1 

^■^'^-0.09 


2.0 X 10" 


8 


- 4.0 X 10" 


5 


Q 1 =7+0.027 


SS 433 


1.53 


0.1 


4 ■sq+O-83 

't-Oy_0.83 


X 10" 


-8 


1 99+0.11 
J^-^^-0.09 


2.0 X 10" 


-8 


- 4.0 X 10" 


-5 


n r,P90.009 


Vela X-1 


5.00 


0.4 




X 10" 


-6 


o.8it^:g^3 


2.0 X 10" 


-8 


- 4.0 X 10" 


-5 


2 556+0155 


GX 301-2 


1.58 


0.8 - 2.0 




X 10" 


-6 


1.86t°J» 


4.0 X 10" 


8 


- 4.0 X 10" 


5 




LMC X-1 


1.36 


~ 23 


9 (-7+0.42 
-1.01 


X 10" 


-7 




5.0 X 10" 


-9 


- 4.0 X 10" 


-5 




SMC X-1 


1.22 


55 


9 9c;+0.9 
•^■^^—0 9 


X 10" 


7 


1 Qq+0.45 


1.0 X 10" 


-8 


- 4.0 X 10" 


'5 


076"^'' °^'' 

'-'•"'"-0.024 


4U 1700-377 


3.83 


~ 0.1 


9 ro+0.43 
^•^°-0.35 


X 10" 


-6 


0.83«:}« 


8.0 X 10" 


-8 


- 4.0 X 10" 


-5 


76+013 



Notes. Lx 

is the X-ray luminosity calculated using respective count rates and source distances listed in Table 1. a is the appropriate power-law 
slope. 
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Fig. 4: The relation bet ween the PDS break frequency and the orbital frequency of the binary system, the findings of 
iGilfanov fc ArefievI (120051 ) regarding short-period LMXRBs being included. The large shaded area toward the bottom of the 
plot is the widest region possible for the dependence fvisc vs. forb according to t he Q-disc, obtained from Eq[T] in the mass ratio 
range 2.91 < g < 22.77 for a = 0.5 and Ra/a - Rudai/a - 0.112 -I- + ^2|f^ (|Gilfanov fc Arefiev|[2005h . The dashed lines are 
predictions for larger values of the disc thickness H/ R with a 
are on the order of marker (or symbol) size. 



1 + 9 ' (1 + >jF 

0.5 and q values indicated on the plot. Most of the errors in ftreak 



3.1. Vela X-1 

Vela X-1 is considered to include a pulsar accreting from 
the violent wind of the supergiant companion. It has been 
observed to exhibit rapid variations in the periodicity of 
its X-ray pulses, which are in dicative of spin-up and spin- 
down phases for its rotation (iDeeter et al.l 119891) and this 
variability is attributed to the torque reversals on the pul- 
sar caused by the change in direction of the gas rotation 
in the temporary accretion disc. This temporary disc has 
been shown, by simulations, to form due to the flip-flop 
instability of the accretion shock fr om the stellar wind to 
the pulsar (Blondin k P ope 200^)- iBovnton et aP (|1984{ ) 
reported a duration of 2-3 days for accretion torque rever- 
sals of Vela X-1, which is close to the viscous time inferred 



from our estimate of break frequency. Hence, the resulting 
disc is very unstable and cannot extend to a large radius 
in the observed 2-3 days as its dispersion is dominated by 
viscosity. 

We estimate the inner and outer radii of a possible ac- 
cretion disc. The magnetic field of the pulsar is rather high, 
^ 10^"^ — 10^^ G, which that means an accretion disc is lim- 
ited b;Y_th£_Alfven_mdius from inside given approximately 
by (jFrank et al.ll2002D 



rAifv 



2.9 X 10^ 



1/7^-2/7^-2/7, ,4/7 



37 M30 



(2) 



where nii = Mn/MQ, i?„ = Rq x 10^ cm, L ^ x 10^^ 
erg/s, and the magnetic moment of the pulsar /i = BnR^ = 
/i30 X 10"^°. Thus, for the Vela pulsar, rAifvm — 3.9 x 10^ cm. 
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Fig. 5: (a)Dependence of the ratio fbreak/ fort on the orbital period of the binary system. Some of the errors in f break / f orb are on 
the order of the marker (or symbol) size. (b)Relation of the normalized ou ter disc radius to the orbital period of the system. The 
arrows represent the possible ranges of the outer disc radii (|Blondinll2000l ). 



There is no distinctive analytical or numerical expression 
for the outer radii of accretion discs of wind-driven HMXBs, 
but the upper limit to a possible disc radius is considered 
to be determined by the tidal forces and Rudai — 0-9Rl 
independent of the characteristics of the disc {Rl is the 
Roche- lobe radius of the compact star). The Roche-lobe 
radii are calcu lated analytically using the expression of 
lEggletonl (Il983h 

Rl ^ 0.49 

a 0.6 + g2/3ln(i + g-i/3)' ^ > 

which implies that Rudai — 6.0 x 10^"'^ cm. In other words, if 
an accretion disc were to form, its radius should be greater 
than the Alfven radius and smaller than the tidal radius. 
This limitation is consistent with our findings since we es- 
timate the outer disc radius to be ~ 4 x 10^ cm using our 
value for the f break / f orb ratio, 2.56, in the theoretical ex- 
pression loT fmsci forb i derived from the Shakura-Sunyaev 
disc model (jGilfanov fc Arefievll2005[) with a thin disc ap- 
proximation. 

We estimated the upper limit to the outer disc radii 
of the other so urces with the same procedure. Simulations 
carried out by iBlondinI ()200Clf ) estimate the disc radius to 
be smaller than our calculations. The limitation depends on 
many parameters such as mass ratio, temperature, Mach 
number, etc. For example, the sizes of warmer discs can be 
smaller by as much as 30%. Nevertheless, these deviations 
would not change the appearance of our results presented 
in the Figl5]-b as a plot of normalized disc radius versus 
orbital period. 

3.2. GX 301-2 

The mass ratio of GX 301-2 is similar to that of Vela X-1, 
which is ^21.7. It has a highly eccentric orbit, e ~ 0.467 
(jSato et al.lll986D with orbital period 41.5 days. The source 
experienced some spin-up episodes since 1984, which are 
interpreted as the signs of the formation of equatorial ac- 
cretion discs from tidal streams found to last for ~ 20 days 



bv lKoh et all (fT997l) . twice in an orbital period. Our anal- 
ysis yielded an inverse break frequency that corresponds to 
~ 11 days. The duration of the temporary accretion disc 
and the timescale we obtained from its break frequency are 
of the same order implying that it has a small accretion 
disc, as Figl5]-b also indicates. 

3.3. 4U 1700-377 

4U 1700-377 is an X-ray eclipsing binary s ystem with an or- 
bital period of 3.4 days (jJones et al.lll973ll . The mass of the 
compact obj ect of the system was found to be ~ 2.44AfQ 
(Clark et al.ll2003) : however, it is still unclear whether it is 
a neutron star or a black hole since no X-ray periodicity has 
been detected that may correspond to the pulse period of 
a neutron star. Nevertheless, it is more common to accept 
that it is a neutron star owing to the similarity of its X-ray 
spectrum to those of accreting neutron stars, and account 
for the lack of pulsations as a result of either a weak mag- 
netic field or an alignment of the magnetic field with the 
spin axis (jClark et al.ll2002[ ). 

The accretion process is considered to be driven by 
a wind flow rather than the Roche-lobe overflow because 
the companion star, HD 153919, under-fills its Roche lobe 
(Conti 1978). Moreover, the ~ 13.8 day variations in the 
X-ray emission and episodic detection of QPOs appear 
to indicate the forma tion of a short-scale transient disc 
()Hong fc HailevI |2004[ ). Its break frequency value corre- 
sponds to ~ 4.4 days. The normalized break frequency, 
0.76, is only a factor of 2-3 smaller then Vela X-1 and GX 
301-2, which indicates that its accretion disc may be of sim- 
ilar size for these sources. This can also be deduced from 
its location on Figl5]-b. 

3.4. SMC X-1 

SMC X-1 is a well-known Roche- lobe-overflow-powered X- 
ray binary ( Papaloizou 197 9) . It is one of the most powerful 
X-ray pulsars and is the only X-ray pulsar displa ying ev- 
idence of steady spin- up (jWoidowski et al.l 119981 ). It also 
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has a super-orbital period of on average ~ 55 days, which 
is assumed to be due t o the precession of a warped disc 
(fWojdowski et aLlll998D . As seen in the figure, its location 
is consistent with other Roche lobe overflow systems. Its 
normalized break frequency is very low at 0.076, implying 
that it has a large accretion disc. 

3.5. OAO 1657-415 

OAO 1657-415 has been proven to be neither a Roche-lobe 
overflow system nor a wind-accreting system although i t 
appears to contain an accretion disc ()Bavkall [19971 l2000h . 
Torque reversals have also been observed to as in the 
case of Vela X-1, but these have timescales of weeks, in 
contrast to the short torque reversal period of Vela X-1 
([Finger fc Prince! 119971 ). The timescale of torque reversal 
and the inverse break frequency ^ 66 days are similar. Its 
normalized break frequency is 0.16, and it should have a 
large disc. On the other hand, that OAO 1657-415 has a 
unique place in the Corbet diagram in-between supergiant 
under-filled Roche Lobe overflow and Be transients suggests 
that there is a longer living accretion disc relative to the 
supergiant wind-fed systems. 

3.6. SS 433 

SS 433 is a black hole binary. It has already been reported 
to have a break in i t s X-r ay spectral band at ~ 10^^ Hz 
bv iRevnivtsev et all ()2006D . Our finding is consistent with 
their break frequency. It is known to have a Roche-lobe- 
overflow-pow ered accretion di sc similar to the low mass X- 
ray binaries (jBegelman et al.ir2006 ').Thus. it is expected to 
be somewhat similar to low mass X-ray binaries with their 
accretion disc features. That is also confirmed by its posi- 
tion on the plot of the correlation of break frequency to or- 
bital frequency (see FigH]). Its normalized break frequency, 
0.05, is close to that of SMC X-1, which indicates that they 
have similar disc properties. 

3.7. LMCX-1 

LMC X-1 is a persistently luminous X-ray source in the 
Large Magellanic Cloud. The estimates for the compact star 
of the binary makes it a good candidate for being a black 
hole binary. It is also unusual because it is a wind driven 
system with a steady state disc. The reason why it is consid- 
ered to be a wind accretor despite having a disc is that the 
optical coun terpart of the com pact star fills only 45% of its 
Roche lobe (jOrosz et al.ll2009D . Since a wind-driven system 
has a disc of smaller radius than a Roche-lobe-overflow bi- 
nary, the disc of LMC X-1 is expected to have a shorter 
viscous timescale than the Roche-lobe overflow systems. 
However, its normalized break frequency, 0.087, is very low 
and its location on the /break- f orb correlation plot is among 
Roche-lobe-overflow binaries despite being expected to be 
near the group of Vela X-1. Ruhlen & Smith (2010) claim 
that LMC X-1 is neither a Roche-lobe-overflow binary nor 
a wind- accretor, but sits at the boundary. Their reasoning 
is that the emission and spectral properties of the source 
do not belong to any specific group. Our result here is con- 
sistent with these findings. A more detailed analysis should 
be performed to account for the unusual behavior of this 
HMXB. 



Figs. [5]-a and [5]-b clearly show the distinction between 
the Roche lobe overflow systems and the wind accreting sys- 
tems. The normalized disc radii presented in FiglS]-b were 
determined from Eqn. ([1]) . In this equation, is taken to 

be the normalized break frequency and is calculated as 
the thickness of standard a-disc 

- 9 V in-2^-l/10n3/20 , ,-21/40 ^ 3/20 pl/8 , 

where Mi. 4 is the mass of the compact object, Mi, in units 
of 1.4Mq, i.e Ml = M1.4 x I.4M0, Rq = lO^Ripkm cm 
is the radius of the compact object, Lx = GMiM / Rq = 
lO^^L^i erg/s is the X-ray luminosity of the source, and 
Rd = lO^^-Rio cm is the outer disc radius. Since the deriva- 
tion of this formula relies on the thin disc approxima- 
tion, our calculations provide lower limits to the disc sizes. 
However, the accretion discs of LMXBs show many signs of 
much thicker accretion discs, eg. by possessing additional 
cor onal flows cove r ing th e thin accretion disc etc., (see 
Gilf anov fc ArefievI (|2005l) for details). The accretion discs 
of HMXBs should have similar structures that make the 
disc seem thicker because of the nature of accretion discs. 
However, we do not expect the irradiation of the disc in high 
mass X-ray binary systems to be as effective as in those of 
low mass X-ray binaries; hence, coronal flow would not ex- 
tend to as large scales as that for H /R ^ 0.1 which the 
case for LMXBs ()Gilfanov fc: ArefievI 120051 ) . Nevertheless, 
the thickening of the accretion disc would not prevent the 
Roche-lobe overflow systems and the wind accreting sys- 
tems being in distinct groups in terms of the low frequency 
break frequencies in their power density spectra. The sub- 
ject of iGilfanov fc ArefievI (|2005^ is the separation of the 
low mass X-ray binaries into two groups. This is not directly 
related to our results although we used a similar method to 
them, because LMXBs are all Roche-lobe overflow systems. 
The grouping they discovered was related to the orbital pe- 
riod of the LMXBs. Their conclusion is that this dichotomy 
is due to the excitation of 3:1 Lindblad resonances in the 
wide systems, which may cause unusual changes in the ac- 
cretion disc, such as the truncation of the disc at a specific 
radius, but does not apply to HMXBs. 

4. Summary 

We have estimated the power density estimates using the 
cosine transform of autocorrelation function, a method that 
can be used for unevenly sampled time series. Using the 
normalized break frequencies of the sources, we have distin- 
guished two groups of high mass X-ray binaries. Those with 
lower normalized break frequencies are more likely to have 
Roche-lobe-overflow accretion discs. These sources are SS 
433, OAO 1657-415, SMC X-1, and LMC X-1. On the other 
hand. Vela X-1, GX 301-2, and 4U 1700-377 have higher 
normalized break frequencies. These sources are wind accre- 
tors. We have discussed the compatibility of the observable 
properties of the sources and normalized break frequencies. 
We have deduced that the value of the normalized break 
frequency is a good estimator for distinguishing wind ac- 
cretors from Roche lobe accretors. This method is especially 
useful for classifying black hole binaries according to their 
accretion mechanisms, since they cannot be shown on the 
Corbet diagram where the spin period of the compact ob- 
ject is compared with the orbital period of the system. 
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